Quantum paramagnetic ground states on the honeycomb lattice and field-induced 

transition to Neel order 
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Bi3Mn40i2(N03) is a recently synthesized spin-3/2 bilayer honeycomb antiferromagnet which 
behaves as a spin liquid down to very low temperatures. Beyond a magnetic field of about 5T, it 
develops long range Neel order. Motivated by this observation, we have studied spin-S* Heisenberg 
models with next neighbor frustrating interactions as well as bilayer couplings on the honeycomb 
lattice. For a model with frustrating second- neighbor exchange, J2, we use a Lindemann-like cri- 
terion within spin wave theory to show that Neel order melts beyond a critical J2. The critical J2 
is found to increase in the presence of a magnetic field, implying the existence of a field-induced 
paramagnet-Neel transition over a range of parameters. For the bilayer model, we use a spin-S* 
generalization of bond operator mean field theory to show that there is a Neel-dimer transition for 
various spin values with increasing bilayer coupling, and that the resulting interlayer dimer state 
undergoes a field induced transition into a state with transverse Neel order. Motivated by a broader 
interest in such paramagnets, we have also studied a spin-3/2 model which interpolates between the 
nearest neighbor Heisenberg model and the Afileck-Kennedy-Lieb-Tasaki (AKLT) parent Hamilto- 
nian. Using exact diagonalization, we have found that there is a single Neel-AKLT quantum phase 
transition in this model. Computing the fidelity susceptibility and assuming a transition in the 0(3) 
universality class, we have located the critical point of this model. In addition, we have obtained 
the spin gap of the AKLT parent Hamiltonian. Our numerics indicate that the AKLT state also 
undergoes a field induced Neel ordering transition. We discuss implications of some of our results 
for experiments on Bi3Mn40i2(N03), and for numerics on the honeycomb lattice Hubbard model. 



The interplay of quantum mechanics and frustrated 
interactions in quantum magnets leads to a variety of 
remarkable phases including spin liquid Mott insulators, 
valence bond crystals, and Bose-Einstein condensates of 
magnonsji Recently, there has been tremendous inter- 
est in novel paramagnetic ground states on the honey- 
comb lattice. A quantum Monte Carlo study of the re- 
pulsive electronic Hubbard model on the honeycomb lat- 
tice has uncovered a spin liquid ground state^ leading to 
a flurry of studies of honeycomb lattice spin liquids. 
Another interesting honeycomb lattice paramagnet is the 
5 = 3/2 AfHeck-Kennedy-Lieb-Tasaki (AKLT) state-^^-^. 
Such an AKLT state is most easily understood by view- 
ing each spin-3/2 as being composed of three spin-1/2 
moments symmetrized on-site, with each spin-1/2 mo- 
ment forming a singlet with one neighbor, leading to a 
ground state which respects all lattice symmetries. This 
state has been suggested as an entanglement resource for 
universal quantum computation!^ Furthermore, an opti- 
cal analogue of the one-dimensional AKLT state has been 
realized in recent experiments,^" raising hopes for alter- 
native realizations of AKLT states in higher dimensions. 

Interest in honeycomb lattice quantum paramagnets 
also stems from experiments^^ on Bi3Mn40i2(N03). The 
crystal field of the MnOg octahedra, together with strong 
Hund's coupling, leads to Heisenberg-like spin-3/2 mo- 
ments on the Mn'^^ ions which form a bilayer honeycomb 
lattice. Despite the bipartite structure, and a large anti- 
ferromagnetic Curie- Weiss constant Qcw ~ —257K, this 
system shows no magnetic order— (or any other phase 



transition) down to T ^ IK. This observation hints 
at frustrating interactions which may lead to interest- 
ing paramagnetic ground statesJ^^— Neutron scattering 
experiments^ on powder samples of Bi3Mn40i2(N03) 
in zero magnetic field indicate that there are short range 
spin correlations in this material, with some antiferro- 
magnetic coupling between the two layers forming the 
bilayer, but negligible interactions between adjacent bi- 
layers. Remarkably, applying a critical magnetic field, 
Be ~ 6 Tesla, leads to sharp Bragg spots consistent with 
three dimensional (3D) Neel orderJ^ 

Motivated by this broad interest in honeycomb lat- 
tice quantum paramagnets, we study various Heisenberg 
models with additional exchange interactions chosen to 
frustrate Neel order. We also consider the effect of a mag- 
netic field on the paramagnetic states which result from 
the destruction of Neel order. We show that applying 
a critical magnetic field to these paramagnetic ground 
states leads to a transition into a state with long range 
Neel order in the plane transverse to the applied field, 
which allows us to make connections with ongoing ex- 
periments and predictions for numerical studies of such 
paramagnetic states. 

We begin with a study of a model with nearest- 
neighbor (Ji) and frustrating second-neighbor (J2) ex- 
change interactions. Such a model is relevant to 
Bi3Mn40i2(N03) as well as numerical studies of the hon- 
eycomb lattice Hubbard model. At the classical level, it 
is known that such frustration leads to the Neel state be- 
coming unstable for J2/J1 > 1/6. In a quantum model 
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with finite S, Neel order is likely to melt for smaller 
J2 / Ji , although the nature of the resulting paramagnetic 
ground state is not known. Sidestepping the issue of what 
state results from quantum melting, we study the mag- 
netic field dependence of the critical J2/J1 required to 
destroy the Neel order. Using spin-wave theory, we show 
that a nonzero magnetic field enhances the critical J2/ Ji, 
opening up a regime where applying a critical field to the 
non-Neel state yields long-range Neel order. 

Next, motivated by the fact that Bi3Mn40i2(N03) 
consists of stacked bilayers, we study a bilayer honey- 
comb magnet where the interlayer exchange interaction 
competes with the intralayer coupling. Using a spin-S* 
generalization22i of the bond operator formalismf2l we 
show that a sufficiently strong bilayer coupling leads to 
an interlayer VBS state. We obtain the Neel to inter- 
layer valence-bond solid (VBS) transition point for var- 
ious spin values, which could be tested using quantum 
Monte Carlo numerics, as well as the triplon dispersion 
in the interlayer VBS. We show that the presence of a 
magnetic field strong enough to overcome the spin gap, 
results in the interlayer VBS undergoing a Bose conden- 
sation transition into a state with long range Neel order 
in the plane transverse to the applied field. For the spin- 
1/2 case, we find that the transition to the interlayer VBS 
state occurs when the interlayer exchange is of the order 
of the interplane exchange (Jc l-3Ji), suggesting that 
spin-1/2 bilayer honeycomb magnets might be a promis- 
ing system to realize this VBS state. For S — 3/2, we 
find that the VBS state is only realized at large interlayer 
couplings, Jc/Ji ?i 6.6. 

Recent attempts to determine the relevant exchange 
couphngs in Bi3Mn40i2(N03) indicate the presence of 
longer range couplin gs in the honeycomb plane J^i^ Fur- 
thermore, while Ref. lis focused on a single layer, the ab 
initio results of Ref. |20| provide evidence for large inter- 
layer couplings. In the light of these reports, our work 
on the J1-J2 model and the bilayer model is perhaps rel- 
evant to the physics of Bi3Mn40i2(N03). Specifically, 
we find that even strong interlayer couplings Jc ^ 2 Ji as 
suggested by Ref. [13 cannot destabilize Neel order in this 
material. We also note that Refs fl9ll20l disagree on the 
sign of the further neighbor couplings within the honey- 
comb layer, suggesting the need for further work on this 
issue. 

Finally, from the viewpoint of broader theoretical in- 
terest, we explore a generalized spin-3/2 model includ- 
ing biquadratic and bicubic spin interactions which in- 
terpolates between a Heisenberg model and the par- 
ent Hamiltonian of the 5 = 3/2 AKLT state. Us- 
ing exact diagonalization, we obtain the spin-gap of the 
AKLT parent Hamiltonian. Wc also compute the fidelity 
susceptibility^^ of this model, and find that it indicates a 
direct AKLT-Neel transition. Using the fidelity suscep- 
tibility and the assumption of an 0(3) critical point, we 
identify the AKLT-Neel transition point in this model. 
By comparing the spin correlations in the singlet ground 
state and in the ground state with 5"*°' = 1, which are 



obtained using the exact diagonalization, we argue that 
a magnetic field applied to the AKLT state results in a 
transition to transverse Neel order. 



I. SECOND-NEIGHBOR EXCHANGE 

It has been suggested that the absence of Neel order 
in Bi3Mn40i2(N03) is linked to non-negligible further 
neighbor interactionsJi We therefore study a minimal 
Hamiltonian, 



J; 



S.-Bj^^f (1) 



where (.) and ((.)) denote nearest and next-nearest neigh- 
bor bonds respectively, and i? is a Zeeman field. Let us 
begin with a classical analysis valid for S = 00. When 
J2 — B ~ 0, the ground state has coUinear Neel order. 
For J2 = and By^O, the spins in the Neel state start 
off in the plane perpendicular to the applied field and 
cant along the field direction until they are fully polar- 
ized for B > 6J1S. For B < 6J1S, the spin components 
transverse to the magnetic field have staggered Neel or- 
der for J2 < Ji/6; for J2 > Ji/6, this gives way to a 
one-parameter family of degenerate (canted) spiralsJ^ 

Incorporating quantum fluctuations is likely to lead to 
melting of Neel order even for J2 < Ji/6. Such fluctua- 
tions are also likely to completely suppress the classical 
spiral order. Using spin wave theory, we argue here that 
a small nonzero B enhances the stability of the Neel or- 
der compared to the zero field case, (i) For small nonzero 
B, spin canting leads to a small decrease, oc i?^, in the 
classical staggered magnetization transverse to the field, 
(ii) On the other hand, one of the two magnon modes 
(labelled f2j^) acquires a nonzero gap oc -B at the F-point 
as shown in Fig. [2] This suppresses low-lying spin wave 
fluctuations. For B<^6JiS, the latter effect overwhelms 
the former, leading to enhanced stability of Neel order. 

Let us discuss this stability line within spin wave the- 
ory. When J2 < Ji/6 and B < 6J1S, Neel order- 
ing is in the plane perpendicular to the magnetic field, 
but the spins also uniformly cant in the direction of 
applied field, to maximally gain Zeeman energy. The 
classical spin state can thus be characterized by Sr = 
S{± cosx, 0, sinx) on the two sublattices. We now define 
new spin operators, denoted by Ti,a, via a sublattice- 
dependent local spin rotation 



,(2) 



where a = 1, 2, is a sublattice index and i sums over each 
unit cell. 

The ground state has all spins pointing toward the new 
local-S*^ axis. To study spin wave fluctuations, we rewrite 
the T operators in terms of Holstein-Primakoff bosons as 
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follows: 



where 



The Hamiltonian can now be rewritten as iJ « Eci + 
Hqu ■ The classical energy Eci is proportional to S"^ , and 
the leading order quantum correction, Hqu, is of order S. 
We get the value of the canting angle x by demanding 
that terms of order S*"^^^, which are linear in the boson 
operators, should vanish, which yields 



smx 



B 
6Ji5' 



The classical energy is given by 



ci 



- Ji cos2x + - J2 



B 
'S 



smx- 



(3) 



(4) 



where N is the number of sites in the honeycomb lat- 
tice. We take the magnetic field B to be of order S, so 
that the Zeeman term ~BSf is treated on the same level 



as the exchange terms JijSi 
correction is given by 



Sj. The leading quantum 



H. 



3NS 



qu 



Ji cos 2x + J2 



NB 



■smx 



k>0 



(5) 



where 



with 



6k, 2 

V &Lk,2 



; H\r 



S X 




(6) 



2J2{coska + cosh + cos{ka + h)} + sinx, 



/k = 3 Ji cos 2x — 6 J2 

ik 

Fk = Ji7kSin2x= |Fk|e*''^ 
Gk = -Ji7kC0S^X: 

where 7k = 1 + c^*'' '' + e-»k (a+fc)^ with unit vectors 
a = x,b = —x/2 + \/3y/2. This Hamiltonian can be di- 
agonalized by a bosonic Bogoliubov transformation. The 
eigenvalues are given by 



= 5VUk±|Fk|)2-|GkP. 



(7) 



The Bogoliubov transformation matrix to rotate into the 
quasiparticle operators is given by 



P = 



C/2X2 
t/2x2 



C'2x2 '5'2x2 
'5'2x2 ^2x2 



u. 



2x2 



1 / — e^'''' e*'"* 

"Til 1 1 



C'2x2— 



cosh 9 
cosh c 



'52x2— 



sinh6' 
sinh ( 



(8) 



, (9) 



where the angles 9 and are given by 

|Gk| 



tanh26' = 
tanh 2(j) ~ 



/k-|Fk|' 

-|Gk| 
/k + li^kl" 



(10) 



The matrix P preserves the commutation relations of the 
bosonic operators and diagonalizes the Hamiltonian, giv- 
ing P^HP = Diag{n-,n+, FigH] shows a plot 
of the magnon dispersion in the Neel state at nonzero B 
along certain high symmetry directions in the Brillouin 
zone. 

The strength of long range magnetic order can be cal- 
culated in this new basis. For example, the in-plane com- 
ponent of the spin is given by 

^ E('^^-i)=(^ + 1/2) cosx- ^ X 

i 

E [cosh26l{ 1-1- 2nB (f^k, -) }+ cosh2(?!){ 1-1- 2nB{r2k,+) }] , ( 1 1 ) 

k>0 

where nBi-) denotes the Bose distribution function. For 
T = 0, we can simply use the 2D Hamiltonian to compute 
this renormalized order parameter. For T ^ 0, we have 
to take into account a small coupling along the third di- 
mension to allow for a stable magnetically ordered state. 
For a layered system with very weak interlayer coupling, 
we can use the 2D Hamiltonian together with an infrared 
cutoff A which is of the order of the interlayer coupling. 
In this case, spin wave modes with energies greater than 
A appear to be 2D spin waves. On the other hand, modes 
with energies below A can be dropped as their contribu- 
tion will be suppressed by phase space factors in the 3D 
problem. In our numerics, we impose this infrared cut- 
off by simply restricting ourselves to a finite system size. 
Finite size automatically cuts off long wavelength modes 
with k < kc ^ 2i: I \In . In our calculations, we have 
restricted our system size to 2 x 120 x 120 spins. This 
corresponds to k^ ^ 0.05, leading to an infrared cutoff of 
A - 0.04JS'. 

As J2 is increased from zero, fluctuations around the 
Neel state increase due to frustration. With increas- 
ing fluctuations, we expect the Neel state to melt when 
fluctuations become comparable to the magnitude of 
the ordered moment. To estimate the 'melting curve', 
we assume that the transverse spin components have 
Neel order along the 5'a;-direction, and use a heuristic 




FIG. 1: Dispersion of magnon modes in the J1-J2 model 
along depicted path in the Brillouin zone for J2=0.15Ji, S— 
3/2 and B^0.5JiS. 
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FIG. 2: r = melting of Neel order for 5=3/2 in the B-J2 
plane (open triangles) obtained using a Lindemann-like crite- 
rion, V (Si) — (Sx)^ = 5{Sx) ■ The region "??" is a quantum 
disordered state - possibly a valence bond solid or a quantum 
spin liquid. Arrow depicts path along which one obtains a 
field-induced transition to Neel order. (Inset) A similar melt- 
ing curve for 5= 1/2. 



Lindemann-like criterion for melting: ^ (5^) — (Sx)^ > 
a{Sx)- The expectation values are evaluated (using lin- 
ear spin wave theory) to order 1, even though the Hamil- 
tonian has terms upto order S only. 

We set a=5 since this leads to melting of Neel order for 
5=1/2 at J2 ~ O.O8J1, in agreement with a recent varia- 
tional Monte Carlo study by Clark et al— The resulting 
Neel melting curves, at zero and nonzero temperatures, 
are shown in Figl2]and Fig|3l 

As shown in Fig l2] and its inset, quantum fluctuations 
at B = lead to melting of Neel order even for J2 < Ji/6 
(i.e., before the classical destruction of Neel order). For 
nonzero B, the 'melting point' moves toward larger J2, 
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FIG. 3: (Color onhne) Melting of Neel order for S= 3/2 in 
the B-J2 plane for depicted nonzero temperatures. To the left 
of the curve, there is stable canted Neel order. To the right, 
combined effects of quantum and thermal fluctuations melt 
the in-plane Neel order. 



leading to a window of J2 over which the quantum disor- 
dered liquid can undergo a field-induced phase transition 
to Neel order. The window of J2 where such physics is 
operative appears to be smah for 5=3/2; however, dis- 
order effects, which tend to suppress the stiffness, may 
enhance this regime. Furthermore, as seen from Fig. |31 
the window of J2 over which we expect field induced Neel 
order is also enhanced at small nonzero temperatures. Fi- 
nally, we expect field induced Neel order even for 5=1/2 
(see inset to Fig. [5]). 

Our results are consistent with recent neutron diffrac- 
tion experiments^^ on Bi3Mn40i2(N03) which find field 
induced Neel order. Our results also explain recent 
Monte Carlo simulations of the classical J1-J2 modeli^ 
with B^O; if J2=0.175 Ji, as in the simulations, increasing 
B at a fixed temperature takes us closer to the melting 
curve as seen from Fig. [3l This may lead to the numeri- 
cally observed enhanced Neel correlations. Nevertheless, 
we expect that there will be no field-induced long-range 
Neel order for J2 = 0.175Ji in the classical model. Fi- 
nally, the Ji— J2 model is a reasonable effective model of 
the insulating phase of the repulsive honeycomb lattice 
Hubbard model, and recent quantum Monte Carlo simu- 
lations find a paramagnetic (spin liquid) insulator over a 
range of repulsion strengths in this model. ^ Our predic- 
tion of field induced Neel order in this paramagnet can be 
verified by including a magnetic field in these quantum 
Monte Carlo simulations. 



II. BILAYER HONEYCOMB LATTICE AND 
THE INTERLAYER DIMER STATE 

The Mn sites in a unit cell of Bi3Mn40i2(N03) form an 
AA stacked bilayer honeycomb lattice. A recent density 
functional theory calculation^" estimates that the inter- 
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layer coupling within each bilayer is large, and may play 
an important role in determining the ground state. If 
the interplane antiferromagnetic exchange Jc is indeed 
large compared to Ji, adjacent spins on the two layers 
could dimerize, leading to loss of Neel order. To study 
this interlayer VBS, we use the Ji— Jc model beginning 
from the limit of Ji = 0; this leads to the spectrum 
E, = -MS{S + 1) - j(j + l)/2), with J = 0, 1, . . . , 25 
denoting the total spin state of the dimer. Restricting 
attention to the low energy Hilbert space spanned by the 
singlet and the triplet states, we define generalized spin-S* 
bond operators via: \s) — s^O), and \a) — tjj|0), where 
|0) is the vacuum, and |q;(= x,y,z)) are related to the 
nij levels of the triplet by \z) — \mj =0), \x) = {\mj = 
-l)-\mj = l))/V2, and |y) =i(|m, = -l) + |m, = 1))/V2. 
Denoting the two spins constituting the dimer, by S^, 
with layer index £ = 0/ 1 , we obtain^^ 



Neel 



VBS 
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Pftlt^, (12) 
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FIG. 4: Phase diagram of the S — 3/2 bilayer honeycomb 
model obtained using bond operator theory, and triplon dis- 
persion along depicted path in the Brillouin zone within the 
interlayer VBS state for Jc/Ji = 7.6 (in units where Ji — l). 



together with the constraint s^s + tl^t^ = I at each site. 

To treat the effect of Ji, we use bond operator mean 
field theory^! which yields a reasonably accurate phase 
diagram for the spin- 1/2 bilayer square lattice Heisen- 
berg modelj ^^'^" Assuming the singlets are condensed in 
the dimer solid, we replace s'^ — s —s, and incorporate 
a Lagrange multiplier in the Hamiltonian which enforces 
(Aj-a) = 1 — on average. Let N be the number of spins 
in each honeycomb layer. We then obtain the Hamilto- 



H 



E 

a,k>0 



-27VC, 



(13) 



describing the dynamics of the triplets. Here 



(4 



ko 



) (with 1, 2 denoting the two sublat- 
tices in each layer) and the matrix Afk takes the form 



r f 



kQl^-kQ2 



/Ak Bk BkN 

BiJ Ak Bl 

Bk Ak Bk 

VS* Bl A^/ 



(14) 



with 



we obtain the mean field values of s and /i which minimize 
the ground state energy subject to the constraint. Solv- 
ing these equations numerically, we find that the spin-S* 
interlayer VBS is a stable phase for Jc > Ji,[S] where 
J43/2]r:!6.6Ji, J^[1]r:!3.5Ji and J^[l/2] « 1.3Ji. Quan- 
tum Monte Carlo studies of this model would be valuable 
in firmly establishing the value of J*[S] as a function of 
S. Fig. 0] shows the triplon dispersion of the S = 3/2 
interlayer VBS state at Jc = 7.6 Ji along high symmetry 
cuts in the hexagonal Brillouin zone. For Jc < J*[S], or 
in the presence of a magnetic field which can close the 
spin gap in the VBS state for Jc > J*[S], the low energy 
triplon mode at the F-point condenses; its eigenvector is 
consistent with Neel order. For the field induced Neel 
state, the Neel ordering is in the plane transverse to the 
applied magnetic field. 

We see Neel order can be destroyed by large bilayer 
coupling, and the resulting interlayer VBS state shows 
field-induced Neel order. However, this is unlikely to be 
the case in Bi3Mn40i2(N03) as it would require a very 
large bilayer coupling Jc ^ 6.6Ji. This ground state 
could be realized in other honeycomb magnetic materials 
with lower spin. 



Ak - Jc-^l-JcS{s + l), (15) 

Bk = '^jkJiS{S + l)s\ (16) 

Here we have defined 7k = 1 + e"**' '' -I- e^*'^ ("+''', with 
unit vectors a = x,b = + \/3y/2, and the constant 

C - - 1 (,s2-l ) - ^{Jc-^i~JcS{S+l)) - i J,s2 5(5+1 ) . (17) 

Diagonalizing this Hamiltonian leads to the ground 
state energy per spin Eg = ^ I]k>o(^k+ +^k-)+C where 
ek± = VA(Ak±2|Bk|). SettmgdEg/ds^ = dEg/dfi = 0, 



III. AKLT VALENCE BOND SOLID 

A particularly interesting spin-gapped ground state of 
a magnet with spin-5 atoms on a lattice of coordina- 
tion number z — 2S, is an AKLT valence bond state. 
Each spin-5 is viewed as being composed of 25 spin- 
1/2 moments symmetrized on-site, with each spin-1/2 
moment forming a singlet with one neighbor;^^^ It was 
originally proposed as an exact realization of Haldane's 
prediction of a spin-gapped ground state in ID inte- 
ger spin systemSf^ Assuming that the Mn^"*" ions in 
Bi3Mn40i2(N03) mainly interact with the three neigh- 
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boring spins in the same plane, this condition is satisfied 
with S = 3/2 and z = 3. The honeycomb lattice AKLT 
state has exponentially decaying spin correlations,— and 
it is the exact, and unique, zero energy ground state of 
the parent Hamiltonian Hakui = Z](ij) Pij*- Here P-^^'' 
denotes a projector on to total spin-^ for a pair of spins 
on nearest neighbor sites (i, j). Denoting Tij = Si ■ Sj, 
we find 



(3) 243 
'J 128 1440 ''^ 



116 ^2 
1440 ''^ 



16 o 
1440 



(18) 



We have investigated, using ED (Exact Diagonaliza- 
tion) on clusters with iV = 12-18 spins, the phase diagram 
of a generalized spin-3/2 model. 



i7Q=(i-g)^s,.s, 



gQH 



AKLT, 



(19) 



which interpolates between a Heisenberg model (at Q=Q) 
and gHAKLT (at Q=l). We set g=1440/243, so that the 
coefficient of S; • Sj in Hq is unity. 

For Q — 0, our analysis of the finite size spectrum 
shows that the ground state energy Eg{N, S^°^), as a 
function of total spin S*^°^, varies as 5'*°'(S'*°* + 1), in 
agreement with the expected Anderson tower for a Neel 
ordered state. It is consistent with earlier wor k^^'^^i^'' 
showing Neel order even for spin- 1/2. To establish the 
Neel- AKLT transition as a function of Q, we study over- 
laps P((5|Q')=l(*s(<3)l*9(Q')>l of the ground state wave 
functions at Q and Q'. As shown in Figl5{a), the overlap 
P{Q\0), of the ground state wavcfunction at Q with the 
Neel state at Q'—O, is nearly unity for Q^O.S, suggesting 
that the ground state in this regime has Neel character. 
For 0.8<(5<1.2, we observe a dramatic drop of P{Q\0) for 
all system sizes, which indicates a Neel-AKLT quantum 
phase transition. 

To locate the transition more precisely, we compute 
the fidelity susceptibilitj^ xf{Q)= 2{1-P{Q\Q+S))/S^, 
with 5^0, which measures the change of the wavcfunc- 
tion when Q — > Q+6. FiglSJb) shows a plot of xf (Q) (with 
(5 = 0.005). We observe a peak in xf{Q) which indicates a 
phase transition; this peak shifts and grows sharper with 
increasing N. Assuming the thermodynamic transition 
is at Q^, and that the peak position Qc{N) satisfies the 
scaling relation {Qc{N) - Q^) N-^/^'' , with i/ « 0.7 
for an 0(3) quantum phase transitio n^^i^^ corresponding 
to triplon condensation, we estimate 0.8. Further 

work is necessary to confirm the nature of the transition. 

The spin gap of Hq, A, (TV) = Eg{N,S'°' = 1) - 
Eg{N,S*-°^ =0), is plotted in FigEfa) for various Q as 
a function of 1/A^. Assuming a finite size scaling form 
AsiN) = Af + b/N, we find a smafi value for for 
Q — 0.0,0.4, consistent with a gapless Neel state, while 
for Q = 0.9, 1.0 we see a robust spin gap as l/N — >• 0. At 
the AKLT point (Q=l), we estimate A^ « 0.6. Given 
our scaling factor g, this yields a value of sa 0.1 for the 
spin gap of ilAKLT- 

Since the spin gap is finite for Q>Qc, we expect that 
applying a critical field Be As will lead to a phase 
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FIG. 5: (Color online) (a) Overlap P{Q\0) of the ground state 
at Q with the Neel state (Q=0) for various system sizes A*', 
showing its rapid drop around the Neel-AKLT transition, (b) 
Fidelity susceptibility Xf{Q) versus Q for various system sizes 
A^, with the peak indicating the Neel-AKLT transition point 
Qc.[N). Inset: Qc (A'') versus 1/iV, together with a fit Qc( A'') = 
Q'^+bN~'^ (with a choice i/ ~ 0.7 assuming an 0(3) quantum 
phase transition in 2D) which leads to Q"^ ^ 0.8. 
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FIG. 6: (Color online) (a) Spin gap, As(Af), versus XjN for 
various Q, with fits to the form A^Af) = AJ° -I- h/N. The 
small values of AJ° for Q = 0.0, 0.4 are consistent with a gap- 
less Neel state. For Q = 0.9, LO, the data are consistent with 
a robust spin gap A^. (b) Si-spin correlations between dis- 
tant sites on the same (^AA) and opposite (AB) sublattices 
for A'^ = 16 system. The spin correlation is Neel-like (±) for 
Q < 1 in the 5*°* = ground state; in the spin gapped AKLT 
state at Q^l, it short ranged and weak in the S*°* = ground 
state but it is strongly enhanced (see arrow) in the lowest 
lying state with Sf* = 1. 



transition; the correlation functions of the lowest lying 
5*°' = 1 state at zero field will then reflect the corre- 
lations of the ground state for > B^- We plot, in 
FiglHKb), the spin correlations on two maximally sepa- 
rated sites (for N = 16) as a function of Q, and make 
the following observations, (i) For S*'"* — 0, the ground 
state also has S"°* = 0, and {S^{i)S^{j)) = (S'^(i)S'^(j)) 
due to spin rotational invariance. At long distance, the 
spin correlation is strong in the Neel phase, but drops 
rapidly to small values upon entering the AKLT state, 
(ii) In the Sl°^ = 1 sector, {S%i)S^ij)) ^ (^,(z)^,(j)). 
Remarkably, in the lowest lying state in this sector, as 
opposed to the 5'*°'=0 ground state, we find a strong en- 
hancement of (only) transverse correlations {Sx{i)Sx{j)) 
between distant sites in the AKLT state; this finite-size 
result suggests that the AKLT state will undergo, beyond 
a critical field, a transition into a state with in-plane Neel 
order. 
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IV. DISCUSSION 

Motivated by recent theoretical and experimental work 
on honeycomb lattice paramagnetic states, we have stud- 
ied various Heisenberg models with competing interac- 
tions. The models we have studied have quantum para- 
magnetic ground states that undergo field-induced phase 
transitions to Neel order. In principle, if such competing 
states are thought to occur in any material, such as in 
Bi3Mn40i2(N03), detailed NMR studies of isolated non- 
magnetic impurities substituted for Mn may help distin- 
guish between these states. For a material with S — 3/2 
moments, the interlayer VBS would have an impurity 
induced 5 = 3/2 local moment on the neighboring site 
in the adjacent layer. The AKLT state would nucleate 
three 5=1/2 moments on neighboring sites in the same 
plane, while spinless impurities in spin gapped Z2 frac- 
tionalized spin liquids, "^"^ do not generically lead to local 
moments. Sharply dispersing triplet excitations expected 
in valence bond solids discussed here could be looked for 
using single-crystal inelastic neutron scattering; by con- 
trast, a spin liquid may not possess such sharp modes. 
Specific heat experiments in a magnetic field could test 
for Bose-condensation of triplet excitations as a route to 
Neel order, which we think describes the transition of the 



AKLT and the interlayer dimer paramagnets. 

If the ground state of Bi3Mn40i2(N03) is a valence 
bond solid, disorder and Dzyaloshinskii-Moriya couplings 
(permitted by the bilayer structure) may be responsible 
for the observed nonzero low temperature susceptibility. 
Finally, dimer crystals with broken symmetry could also 
be candidate ground states in Bi3Mn40i2(N03); if this 
is the case, disorder must be responsible for wiping out 
the thermal transition expected of such crystals. These 
are interesting directions for future research. 
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